Mechanistic modeling for optimal design of dissolvable microneedle-based patches for transdermal drug delivery.

Polymeric microneedle (MN)-based patches are an efficient, non-invasive, and painless means of drug delivery through the skin to systemic circulation. The design of these MN-based patches can be customized for various drug delivery applications, particularly modified release of drugs. In this study, we developed a mathematical model of drug delivery via MN-based patches to study the effect of patch design properties on drug delivery kinetics and systemic pharmacokinetics (PK). We calibrated the model against two representative formulations: a rapid release patch of naloxone and a sustained-release patch of levonorgestrel. The model was then applied to assess the relative significance of model parameters in governing systemic PK of drugs and obtain optimal design parameters to achieve therapeutically meaningful drug levels in a clinical setting. We identified the importance of drug loading fraction, MN base radius, and MN height as the key control parameters responsible for drug PK.Clinical Relevance- Through the application of modeling and simulation, we can improve drug delivery from MN-based patches by identifying optimal design parameters to support the clinical translation of these novel drug delivery systems.

Hybrid Cellular Automata Modeling Reveals the Effects of Glucose Gradients on Tumour Spheroid Growth.

PURPOSE: In recent years, mathematical models have become instrumental in cancer research, offering insights into tumor growth dynamics, and guiding the development of pharmacological strategies. These models, encompassing diverse biological and physical processes, are increasingly used in clinical settings, showing remarkable predictive precision for individual patient outcomes and therapeutic responses. METHODS: Motivated by these advancements, our study introduces an innovative in silico model for simulating tumor growth and invasiveness. The automated hybrid cell emulates critical tumor cell characteristics, including rapid proliferation, heightened motility, reduced cell adhesion, and increased responsiveness to chemotactic signals. This model explores the potential evolution of 3D tumor spheroids by manipulating biological parameters and microenvironment factors, focusing on nutrient availability. RESULTS: Our comprehensive global and local sensitivity analysis reveals that tumor growth primarily depends on cell duplication speed and cell-to-cell adhesion, rather than external chemical gradients. Conversely, tumor invasiveness is predominantly driven by chemotaxis. These insights illuminate tumor development mechanisms, providing vital guidance for effective strategies against tumor progression. Our proposed model is a valuable tool for advancing cancer biology research and exploring potential therapeutic interventions.

Mechanistic modeling-guided optimization of microneedle-based skin patch for rapid transdermal delivery of naloxone for opioid overdose treatment.

Naloxone, an FDA-approved opioid inhibitor, used to reverse opioid overdose complications has up till date faced challenges associated with its delivery. Limitations include the use of invasive delivery forms and the need for frequent redosing due to its short half-life. The goal of the current study was to design a transdermal rapidly dissolving polymeric microneedle (MN) patch with delivery and pharmacokinetic properties comparable to that seen with the commercially available NAL products, eliminating their delivery limitations. Patches of varying dimensions (500 µm; 100 array,800 µm; 100array, and 600 µm; 225 array) were fabricated to evaluate the effect of increasing MN length, and density (no. of needles/unit area) on drug release. Drug dose in each of these patches was 17.89 ± 0.23 mg, 17.2 ± 0.77 mg, and 17.8 ± 1.01 mg, respectively. Furthermore, the insertion efficiency of each of the MN patches was 94 ± 4.8%, 90.6 ± 1.69%, and 96 ± 1.29%, respectively. Compared to passive permeation, a reduced lag time of about 5-15 min was observed with a significant drug flux of 15.09 ± 7.68 g[Formula: see text]/cm2/h seen in the first 1 h (p < 0.05) with the array of 100 needles (500 µm long). Over 24 h, a four and ten-fold increase in permeation was seen with the longer length and larger density MN patch, respectively, when compared to the 500 µm (100 array) patch. Model simulations and analyses revealed the significance of needle base diameter and needle count in improving systemic pharmacokinetics of NAL.

Translational Modeling Identifies Synergy between Nanoparticle-Delivered miRNA-22 and Standard-of-Care Drugs in Triple-Negative Breast Cancer.

PURPOSE: Downregulation of miRNA-22 in triple-negative breast cancer (TNBC) is associated with upregulation of eukaryotic elongation 2 factor kinase (eEF2K) protein, which regulates tumor growth, chemoresistance, and tumor immunosurveillance. Moreover, exogenous administration of miRNA-22, loaded in nanoparticles to prevent degradation and improve tumor delivery (termed miRNA-22 nanotherapy), to suppress eEF2K production has shown potential as an investigational therapeutic agent in vivo. METHODS: To evaluate the translational potential of miRNA-22 nanotherapy, we developed a multiscale mechanistic model, calibrated to published in vivo data and extrapolated to the human scale, to describe and quantify the pharmacokinetics and pharmacodynamics of miRNA-22 in virtual patient populations. RESULTS: Our analysis revealed the dose-response relationship, suggested optimal treatment frequency for miRNA-22 nanotherapy, and highlighted key determinants of therapy response, from which combination with immune checkpoint inhibitors was identified as a candidate strategy for improving treatment outcomes. More importantly, drug synergy was identified between miRNA-22 and standard-of-care drugs against TNBC, providing a basis for rational therapeutic combinations for improved response CONCLUSIONS: The present study highlights the translational potential of miRNA-22 nanotherapy for TNBC in combination with standard-of-care drugs.

A Multiscale Model to Identify Limiting Factors in Nanoparticle-Based miRNA Delivery for Tumor Inhibition.

MicroRNA-based gene therapy for cancer treatment via nanoparticles (NPs) requires navigation of multiple physical and physiological barriers in order to efficiently deliver the miRNAs to the cancer cell cytoplasm. We here present a mathematical model to investigate the variability associated with tumor, NP, and miRNA characteristics, and identify the limiting factors in miRNA delivery to tumors. Through global parameter analysis, the miRNA release rate from NPs and NP degradability were found to have the most significant impact on cytosolic accumulation of miRNAs. These NP properties can be fine-tuned in order to optimize the delivery system for enhancing the efficacy of miRNA-based therapy.Clinical Relevance-Understanding the effect of nanoparticle, tumor, and miRNA characteristics in governing the efficacy of miRNA-based cancer therapy will support its clinical translation.

Innate Immunity Plays a Key Role in Controlling Viral Load in COVID-19: Mechanistic Insights from a Whole-Body Infection Dynamics Model.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a pathogen of immense public health concern. Efforts to control the disease have only proven mildly successful, and the disease will likely continue to cause excessive fatalities until effective preventative measures (such as a vaccine) are developed. To develop disease management strategies, a better understanding of SARS-CoV-2 pathogenesis and population susceptibility to infection are needed. To this end, mathematical modeling can provide a robust in silico tool to understand COVID-19 pathophysiology and the in vivo dynamics of SARS-CoV-2. Guided by ACE2-tropism (ACE2 receptor dependency for infection) of the virus and by incorporating cellular-scale viral dynamics and innate and adaptive immune responses, we have developed a multiscale mechanistic model for simulating the time-dependent evolution of viral load distribution in susceptible organs of the body (respiratory tract, gut, liver, spleen, heart, kidneys, and brain). Following parameter quantification with in vivo and clinical data, we used the model to simulate viral load progression in a virtual patient with varying degrees of compromised immune status. Further, we ranked model parameters through sensitivity analysis for their significance in governing clearance of viral load to understand the effects of physiological factors and underlying conditions on viral load dynamics. Antiviral drug therapy, interferon therapy, and their combination were simulated to study the effects on viral load kinetics of SARS-CoV-2. The model revealed the dominant role of innate immunity (specifically interferons and resident macrophages) in controlling viral load, and the importance of timing when initiating therapy after infection.

Imaging-Based Subtypes of Pancreatic Ductal Adenocarcinoma Exhibit Differential Growth and Metabolic Patterns in the Pre-Diagnostic Period: Implications for Early Detection.

BACKGROUND: Previously, we characterized subtypes of pancreatic ductal adenocarcinoma (PDAC) on computed-tomography (CT) scans, whereby conspicuous (high delta) PDAC tumors are more likely to have aggressive biology and poorer clinical outcomes compared to inconspicuous (low delta) tumors. Here, we hypothesized that these imaging-based subtypes would exhibit different growth-rates and distinctive metabolic effects in the period prior to PDAC diagnosis. MATERIALS AND METHODS: Retrospectively, we evaluated 55 patients who developed PDAC as a second primary cancer and underwent serial pre-diagnostic (T0) and diagnostic (T1) CT-scans. We scored the PDAC tumors into high and low delta on T1 and, serially, obtained the biaxial measurements of the pancreatic lesions (T0-T1). We used the Gompertz-function to model the growth-kinetics and estimate the tumor growth-rate constant (α) which was used for tumor binary classification, followed by cross-validation of the classifier accuracy. We used maximum-likelihood estimation to estimate initiation-time from a single cell (10-6 mm3) to a 10 mm3 tumor mass. Finally, we serially quantified the subcutaneous-abdominal-fat (SAF), visceral-abdominal-fat (VAF), and muscles volumes (cm3) on CT-scans, and recorded the change in blood glucose (BG) levels. T-test, likelihood-ratio, Cox proportional-hazards, and Kaplan-Meier were used for statistical analysis and p-value <0.05 was considered significant. RESULTS: Compared to high delta tumors, low delta tumors had significantly slower average growth-rate constants (0.024 month-1 vs. 0.088 month-1, p<0.0001) and longer average initiation-times (14 years vs. 5 years, p<0.0001). α demonstrated high accuracy (area under the curve (AUC)=0.85) in classifying the tumors into high and low delta, with an optimal cut-off of 0.034 month-1. Leave-one-out-cross-validation showed 80% accuracy in predicting the delta-class (AUC=0.84). High delta tumors exhibited accelerated SAF, VAF, and muscle wasting (p <0.001), and BG disturbance (p<0.01) compared to low delta tumors. Patients with low delta tumors had better PDAC-specific progression-free survival (log-rank, p<0.0001), earlier stage tumors (p=0.005), and higher likelihood to receive resection after PDAC diagnosis (p=0.008), compared to those with high delta tumors. CONCLUSION: Imaging-based subtypes of PDAC exhibit distinct growth, metabolic, and clinical profiles during the pre-diagnostic period. Our results suggest that heterogeneous disease biology may be an important consideration in early detection strategies for PDAC.

Innate immunity plays a key role in controlling viral load in COVID-19: mechanistic insights from a whole-body infection dynamics model.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a pathogen of immense public health concern. Efforts to control the disease have only proven mildly successful, and the disease will likely continue to cause excessive fatalities until effective preventative measures (such as a vaccine) are developed. To develop disease management strategies, a better understanding of SARS-CoV-2 pathogenesis and population susceptibility to infection are needed. To this end, physiologically-relevant mathematical modeling can provide a robust in silico tool to understand COVID-19 pathophysiology and the in vivo dynamics of SARS-CoV-2. Guided by ACE2-tropism (ACE2 receptor dependency for infection) of the virus, and by incorporating cellular-scale viral dynamics and innate and adaptive immune responses, we have developed a multiscale mechanistic model for simulating the time-dependent evolution of viral load distribution in susceptible organs of the body (respiratory tract, gut, liver, spleen, heart, kidneys, and brain). Following calibration with in vivo and clinical data, we used the model to simulate viral load progression in a virtual patient with varying degrees of compromised immune status. Further, we conducted global sensitivity analysis of model parameters and ranked them for their significance in governing clearance of viral load to understand the effects of physiological factors and underlying conditions on viral load dynamics. Antiviral drug therapy, interferon therapy, and their combination was simulated to study the effects on viral load kinetics of SARS-CoV-2. The model revealed the dominant role of innate immunity (specifically interferons and resident macrophages) in controlling viral load, and the importance of timing when initiating therapy following infection.

Image-guided mathematical modeling for pharmacological evaluation of nanomaterials and monoclonal antibodies.

While plasma concentration kinetics has traditionally been the predictor of drug pharmacological effects, it can occasionally fail to represent kinetics at the site of action, particularly for solid tumors. This is especially true in the case of delivery of therapeutic macromolecules (drug-loaded nanomaterials or monoclonal antibodies), which can experience challenges to effective delivery due to particle size-dependent diffusion barriers at the target site. As a result, disparity between therapeutic plasma kinetics and kinetics at the site of action may exist, highlighting the importance of target site concentration kinetics in determining the pharmacodynamic effects of macromolecular therapeutic agents. Assessment of concentration kinetics at the target site has been facilitated by non-invasive in vivo imaging modalities. This allows for visualization and quantification of the whole-body disposition behavior of therapeutics that is essential for a comprehensive understanding of their pharmacokinetics and pharmacodynamics. Quantitative non-invasive imaging can also help guide the development and parameterization of mathematical models for descriptive and predictive purposes. Here, we present a review of the application of state-of-the-art imaging modalities for quantitative pharmacological evaluation of therapeutic nanoparticles and monoclonal antibodies, with a focus on their integration with mathematical models, and identify challenges and opportunities. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease Diagnostic Tools > in vivo Nanodiagnostics and Imaging Nanotechnology Approaches to Biology > Nanoscale Systems in Biology.

Mathematical Modeling to Address Challenges in Pancreatic Cancer.

Pancreatic Ductal Adenocarcinoma (PDAC) is regarded as one of the most lethal cancer types for its challenges associated with early diagnosis and resistance to standard chemotherapeutic agents, thereby leading to a poor five-year survival rate. The complexity of the disease calls for a multidisciplinary approach to better manage the disease and improve the status quo in PDAC diagnosis, prognosis, and treatment. To this end, the application of quantitative tools can help improve the understanding of disease mechanisms, develop biomarkers for early diagnosis, and design patient-specific treatment strategies to improve therapeutic outcomes. However, such approaches have only been minimally applied towards the investigation of PDAC, and we review the current status of mathematical modeling works in this field.